Chandra Observations of Six QSOs at z« 3 
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ABSTRACT 

We report the results of our Chandra observations of six QSOs at z ~ 3 from the Palomer 
Transit Grism Survey. Our primary goal is to investigate the possible systematic change of a ox 
between z > 4 and z ~ 3, between which a rapid rise of luminous QSO number density with 
cosmic time is observed. The summed spectrum showed a power-law spectrum with photon index 
of r « 1.9, which is similar to other unabsorbed AGNs. Combining our z ~ 3 QSOs with X-ray 
observations of QSOs at z > 4 from literaure/archive, we find a correlation of a x with optical 
luminosity. This is consistent with the fact that the luminosity function slope of the luminous 
end of the X-ray selected QSOs is steeper than that of optically-selected QSOs. We discuss an 
upper limit to the redshift dependence of a ox using a Monte-Carlo simulation. Within the current 
statistical errors including the derived limits on the redshift dependence of a ox , we found that 
the behaviors of the X-ray and optically-selected QSO number densities are consistent with each 
other. 



Subject headings: galaxies: active — galaxies: high-redshift 
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1. Introduction 

Large-scale optical surveys show that the lumi- 
nous QSO number density peaks at 1.7 < z < 2.7, 
before which (in cosmic time) the QSO number 
density grows rapidly and after which the den- 
sity steadily decays until the present epoch (e.g. 
Boyle et al. 1988; Warren, Hewett, & Osmer 1994; 
Schmidt, Schneider, & Gunn 1995). The first 
indication of a density decline (with redshift) at 



1 Present Address:LPS-Berlin, Kopernickerstr. 325, 
12555, Berlin, Germany 
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z > 3 was reported in the pioneering work of Os- 
mer (1982). The Palomar Transit Grism Survey 
(PTGS) (Schneider, Schmidt, & Gunn 1994) was 
designed to investigate this possible "redshift cut- 
off', and produced a sample of 90 z > 2.7 QSOs; 
an analysis of the PTGS Schmidt, Schneider, & 
Gunn (1995) (hereafter SSG95) revealed a very 
rapid growth, by a factor of ~ 3 — 5, of the number 
density of luminous QSOs from z ~ 4.5 to ~ 2.7. 
This result was verified by Fan et al. (2001) us- 
ing a sample of 3.6 < z < 5.0 quasars found in 
the Sloan Digital Sky Survey (SDSS; York et al. 
2000) . This redshift cutoff has also been found in 
the number densities of radio flat-spectrum QSOs 
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(Shaver et al. 1996; Wall et al. 2005). 

In the X-ray regime, however, the size of a sam- 
ple based on the combination of various ROSAT 
surveys, Miyaji, Hasinger, & Schmidt (2000) (her- 
aaftcr, SXLF1) was not sufficient to accurately 
probe this "growth" phase of the number density 
of X-ray-selected luminous (log L x > 44.5) QSOs. 
While this study appears to show a flat number 
density in z > 2.7, the uncertainties produced by 
the small number of quasars limited the robustness 
of the conclusions. More recent studies including 
the results from Chandra and XMM-Newton sur- 
veys (Cowie et al. 2003; Fiore et al. 2003; Ueda et 
al. 2003; Barger et al. 2005) showed that the num- 
ber density curve is luminosity-dependent and a 
trend that as luminosity goes lower, the density 
peak shifts to lower redshifts (as some authors 
call an anti-hierarchical AGN evolution and others 
call a down-sizing of the AGN activity). However, 
these studies did not sample well the high-redshift, 
high-luminosity regime, in which optically selected 
QSOs show a decline at z > 3. 

In a recent study based on an updated soft X- 
ray sample including Chandra and XMM-Newton 
surveys, Hasinger, Miyajife Schmidt (2005) (here- 
after SXLF3) also measured the number densities 
in various luminosity bins with a better accuracy 
over a large range in the redshift-luminosity space. 
The study revealed this early "growth" of AGN 
number density (or decline with increasing z) at 
z> 3 for AGN/QSO luminosities at (log L x < 44- 
45). The results from "Champs" survey, which 
was designed to optimally trace the high-redshift, 
high-luminosity (logL x > 44.5) regime with im- 
proved statistics, found that the density curves X- 
ray selected QSOs (\ogL x > 44.5) declines with 
increasing redshift at z > 3. This decline is, how- 
ever, shallower than that seen in optically-selected 
QSOs. 

The dependence of the optical to X-ray flux ra- 
tio (customarily expressed by the quantity cn ox , 
the effective spectral index between the rest-frame 
2500 A and 2 keV) 2 on redshift and luminosity 
has been a key issue in X-ray observations of high- 
redshift QSOs and has important implications for 
possible differences in the AGN evolution traced 
by X-ray and optical samples. Some authors (e.g. 
Vignali, Brandt, & Schneider 2003b; Strateva et 

2 We use /„ oc u axx , where XX is any subscript to a. 



al. 2005) found that a ox strongly depends on lu- 
minosity, with Z/ x oc £/ojytf ' and no evidence for 
any evolution of the X-ray properties with red- 
shift. Bcchtold et al. (2003) found that variations 
in a ox depends primarily on redshift. The depen- 
dence may be sensitive to the selection effects, in- 
cluding but not limited to whether the sample is 
optically-selected or X-ray selected. Yuan et al. 
(1998b) pointed out that such an apparent non- 
linearity of the luminosity correlation in two bands 
can arise from the difference in the luminosity vari- 
ations in the two bands. 

In order to investigate the redshift dependence 
of the optical-to-X-ray luminosity ratios and its 
impact on the density curves of luminous QSOs 
in X-ray and optically selected samples at z > 
3, we have obtained Chandra observations of six 
PTGS QSOs with redshifts between 2.91 and 2.96. 
This is the era of maximum number density of 
luminous QSOs. There was practically no sys- 
tematic observations in X-rays in this redshift 
regime before. Thus our observations also serve 
to fill this observation gap. Throughout this 
paper we adopt (_ff [kms _1 Mpc -1 ], fi m , ft a) = 
(70/i7o,0.3,0.7) and /170 = 1 unless otherwise 
noted. 

2. Observations and Analysis 
2.1. The Sample and Observation 

The original motivation of the program was to 
compare the mean a ox values of z ~ 3 and z > 4 
QSOs. In Chandra Cycle 4, six z ~ 3 QSOs from 
the PTGS (out of 15 QSOs proposed) have been 
observed. None of these QSOs are broad absorp- 
tion line (BAL) QSOs. Table 1 shows the observed 
targets, log of observations, optical AB magni- 
tudes at the object's rest frame of 1450 A. The core 
radio loudness R L = log 10 f v [20cm] //„ [4400 A] is 
also shown, where those with i?L < 1 and > 1 are 
divided into radio-quiet QSOs (RRQ) and radio- 
loud QSOs (RLQ) respectively (Wilkes 2000 and 
references therein). The fluxes are in the object's 
rest frame, calculated assuming radio and optical 
spectral indices of a r = —0.8 and a a — —0.79 
respectively. The radio data are from the NVSS 
(PC 0041+0024 Condon et al. 1998) or FIRST (all 
others Becker, White, & Helfand 1995) surveys. 
Only one QSO (PC 1035+4747) was detected in 
the radio band and for others, 3tr upper limits of 
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i?L are shown. The only radio-detected QSO, PC 
1035+4747, has R h = 2.3; this object falls well 
into the RLQ regime. Two of the Rh upper limits 
(1.5 for PC 0041+0214 and 1.2 for PC1000+4751) 
are above the RQQ/RLQ border, but their limits 
are well below the peak of the Rh distribution of 
RLQ. Thus we tentatively classify them as RQQs. 

All observations have been made with the 
Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003) and the targets aimed at 
the default position of ACIS-S with the S3-chip. 
In all cases except PC 0947+5628, X-ray counter- 
parts have been found within 1.5" of the cataloged 
optical center of the target QSOs, consistent with 
the combined systematic error on the absolute as- 
trometries of l"-2" in both the Chandra data prod- 
ucts and the PTGS survey. Although the detected 
X-ray source closest to PC0947+5628 was 2.8" 
away from the catalogued position of the QSO, 6 
other X-ray sources in the same ACIS observation 
also had optical counterparts at ~ 2.8" away with 
practically the same offset directions. Thus we 
also identify the X-ray source with PC0947+5628. 

2.2. Individual and Stacked Spectra 

For all six observations, we have extracted the 
pulse-height spectra using an extraction radius of 
2". Spectra, response matrices (rmf) and ancillary 
response files (arf ) were created using the software 
package CIAO 3.0.2 or later versions, in conjunc- 
tion with the calibration database CALDB 2.26 
or later versions. These versions enabled construc- 
tion of the response files which takes the time- 
dependent low energy efficiency degradation into 
account. Due to low number counts of the in- 
volved objects, changes due to further updates of 
the calibration have negligible effects. The spec- 
tral analyses were made to the pulse-height chan- 
nels corresponding to observed photon energies of 
0.3-7 keV. Background level is typically ~ 0.05 
counts in the extraction radius and is thus neg- 
ligible. The spectral analysis were made with 
XSPEC 11.2. In spite of small number of X- 
ray photons, the negligible background and the 
use of the XSPEC implementation of the Cash 
(1979) C-statistics allowed placement of some con- 
straints on the spectral indexes, although there are 
not sufficient number of photons in any individual 
spectrum to simultaneously constrain the intrin- 
sic absorption column density. The results of the 



spectral fits with a single power-law with a photon 
index T and the Galactic absorption N§ 20 [cm -2 ] 
Dickey & Lockman (1990) at the position of the 
QSOs are shown in Table 2. We see that the only 
radio-detected QSOs photon index of T = 0.4 + 1.0 
is constrained to be harder than the mean QSO 
spectrum. The rest-frame 2-10 keV luminosities 
(log L x ; logarithm is base- 10) are also shown as 
well as the source counts in 0.3-7 keV. The rest- 
frame 2-10 keV corresponds to observed frame 0.5- 
2.5 keV. The B-band absolute magnitudes (Mb), 
that have been recalculated using our default cos- 
mology (Sect. 1) and the optical spectral index 
a Q = -0.79 (following Vignali et al. 2003c; Fan et 
al. 2001), are also listed here. Using a = —0.5 
(Schneider et al. 2001) increases Mb by 0.35. 




channel energy (keV) 



Fig. 1. — The stacked pulse-height spectrum of the 
6 QSOs at z ~ 3 with the folded best-fit power- 
law model and residuals in terms of the data-to- 
model ratio. The spectrum is rebinned for display, 
but the actual fit was made to a higher resolution. 
Confidence contours at levels of AC = 2.3 and 
4.6, 9.2 (corresponding to the 68%, 90% and 99% 
confidence levels for two interesting parameters re- 
spectively) in the Ah22 — T plane are also shown. 
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We also analyzed the summed spectrum from 
all six QSOs. Because of the small spread of the 
redsifts of these QSOs, we can analyze the summed 
spectrum assuming a single redshift. The response 
matrix for the summed spectrum was constructed 
by a source-count weighted mean the 6 matrices. 
The Galactic column densities N^ 20 of these six 
QSOs are 0.9, 1.0, 1.0, 1.3, 2.0 and 2.8. The 
fit was made to a model with the sum of three 
power-laws, with different Galactic absorptions of 
^H20 = 1-0,2.0 and 2.8 respectively, where the 
four QSOs with 0.9 < N§ 20 < 1.3 were repre- 
sented by a single column density of 1.0. The 
photon indices of all the three components were 
set to equal and the ratios of the 3 normalizations 
were fixed to those of the total source counts of 
the QSOs with N§ 20 of 0.9-1.3, 2.0 and 2.8 re- 
spectively. Also an intrinsic absorption compo- 
nent Nft 2 Q is included with z = 2.93, which is a 
source-count weighted mean redshift of the sam- 
ple. Again, the C-statistics was used for the fit. 
The summed spectrum is well represented by a 
single power-law with T ~ 1.9 and no intrinsic ab- 
sorption, as shown in the last entry of Table 2. 
The pulse-height spectrum, folded best-fit power- 
law model, and fit residuals in terms of data-to- 
model ratio are shown in Fig. 1, with confidence 
contours for the intrinsic absorption versus pho- 
ton index space. Removing the one RLQ (PC 
1035+4747) from the analysis did not change the 
fitted parameters and errors to the smallest dig- 
its displayed in Table 2. This result is consistent 
with the mean slopes of RQQs and unabsorbed 
AGNs measured in the 2-10 keV in the rest frame 
over a wide range of redshift and luminosity (e.g. 
Vignali et al. 2003a, c). Note, however, that the 
stacked spectrum is dominated by a few brightest 
sources, with 60% of the photons coming from the 
two brightest objects. The quoted error only in- 
cludes the statistical error of photon counts. The 
sampling error is estimated by a bootstrapping 
method, where the 90% error range was deter- 
mined by 500 bootstrap runs of a photon-count 
weighted mean best-fit T values. The results were 
(90% bootstrap errors) (r) = 1.85j£|i (1.97±0.29 
with the RLQ removed). 

2.3. The Optical to X-ray Index (a ox ) 

The optical (rest-frame ultraviolet) to X-ray 
flux ratio of a QSO is customarily expressed in 



terms of the effective index a ox between 2 keV 
and 2500A in the QSOs rest frame. Upon calcu- 
lating a ox , we assumed T = 2.0 for all, which is 
the average QSO spectral index. The result of the 
spectral analysis of all but one is consistent with 
the canonical spectral index of T = 1.9-2.0. The 
RLQ PC 1035+4747 has T ~ 0.6 ± 0.8 and using 
F = 0.6 for K-correction gives a ox = —2.06. Other 
than this one, the main source of errors in a ox is 
the X-ray flux. Even for the source with the small- 
est source count (PC 1035+4747) the la error on 
a ox is ~ 0.05. A decrease of F by 0.2 leads to an 
decrease of a ox by 0.04 at z <~ 3. Using a Q = —0.5 
instead of —0.79, a ox increases by 0.03. 

3. Redshift and Mb Dependences of a ox 

This program is mainly focused on the system- 
atic difference in a ox between z ~ 3 and z > 4. 
Because a ox values of QSOs show a large scatter 
and we only sample a small number of QSOs in 
the redshift-luminosity regimes of our interest, the 
sampling error is the dominant effect in the error 
budget of the mean value ({a ox )) of QSOs, which 
can be estimated by a/ a/-/Vq, where a is the stan- 
dard deviation of the a ox distribution of Nq QSOs. 
Note that the a/ \/Nq estimation of the standard 
deviation of the mean is also valid for small Nq . It 
is well known that this estimator gives the exact 
confidence range of Gaussian la when the parent 
a ox distribution is a Gaussian and it is widely used 
in more general cases. The a ox distributions in 
Yuan et al. (e.g. 1998a); Strateva et al. (e.g. 2005); 
Vignali, Brandt, & Schneider (e.g. 2003b) are well 
characterized by a Gaussian, which justifies the 
use of this estimator in our analysis. For our sam- 
ple, we obtain (a ox ) = -1.65 + 0.05 (-1.62 + 0.05, 
for the 5 RQQs only). Our results are compared 
with those of z > 4 QSOs in Table 4. Ideally, 
we would like to compare with z > 4 QSOs in 
the same luminosity range (-27. 0> Mb > -26.3) 3 . 
Unfortunately all but a few of the z > 4 QSOs ob- 
served previously with X-rays found in literature 
(Bechtold et al. 2003; Vignali et al. 2001, 2003a, 
2005) are more luminous (Mb < —27) than those 
in our sample. Keeping this limitation in mind, 

3 The magnitude limit used by SSG95 of M B = -26.0, who 
used h7o = 5/7, S7 m = 1, S^a = with a a = —0.5, corre- 
sponds to Mb = —26.47 and -26.55 for our default cosmol- 
ogy and a = —0.79 at z = 3 and 5 respectively. 
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we compare our result with the mean a ox val- 
ues in Tables 3 and Al in Vignali et al. (2003c) 
and Table 3 in Vignali et al. (2005), who used 
the same methods and assumptions as our work 
in deriving fluxes and a ox values. Their Table Al 
includes recalculated a ox values of z > 4 QSOs 
previously observed by Chandra (including those 
in Bechtold et al. 2003) using the same method. 
Their recalculations also took the time-dependent 
low-energy degradation of ACIS-S quantum effi- 
ciency into account. BALQSOs SDSS 1129-0142 
and SDSS 1605-0122 (Vignali et al. 2003a) have 
been excluded from the following analysis, because 
BALQSO are known to be X-ray weak due to a 
heavy absorption. Furthermore, we have derived 
X-ray fluxes of further three z > 4 QSOs from pub- 
lic archival Chandra ACIS-S data and calculated 
their a ox values in the same way, as summarized 
in Table 3. The 1450 Amagnitudes of these three 
have been obtained from SSG95 (PC 0910+5625) 
or the spectra from the SDSS DR3 database 4 af- 
ter corrections for Galactic extinction (the others) . 
These three have been included in our statistical 
analysis. 

The combined z > 4 sample is somewhat het- 
erogeneous, as it consists of data obtained by a 
variety of programs, each with different interests 
and strategies. In our statistical analysis, we have 
divided the sample into five groups (Group A-E, 
with Group A being our z <~ 3 sample) as shown 
in Table 4. Figure 2(b) shows the scatter diagram 
of the combined sample in the z-M-q plane with 
symbols showing the group membership. For each 
group, we have calculated the mean value (a ox ) 
and the standard deviation of the mean. The 
combined sample includes six upper limit a ox val- 
ues (no X-ray detection) and we have used the 
Kaplan-Meier estimator to calculate the mean and 
standard deviation of the mean using the ASURV 
software(Feigelson & Nelson 1985), when neces- 
sary. Figure 2(a) and (c) plots (a ox ) of each group 
against z and Mb respectively. 

Figure 2(b) shows a rather incoherent scatter, 
reflecting the interests of individual original ob- 
serving programs. As a result, z and Mb do not 
show strong monotonic z versus luminosity cor- 
relation, unlike flux-limited samples. This is for- 
tunate, because if such a correlations existed, it 

4 http: / /www. sdss.org/ 



would be difficult to separate he redshift and lu- 
minosity dependences of a ox . 

Figs. 2(a) and (c) seems to support the conclu- 
sion by (Vignali, Brandt, & Schneider 2003b) that 
the a ox is anti-correlated with ultraviolet luminos- 
ity, while the correlation with z is weak, if any. We 
made a linear regression analysis of the dependent 
parameter a ox against the two independent pa- 
rameters, namely z and Mb- Wc also made sepa- 
rate one-independent parameter regressions for z 
or Mb versus the dependent parameter a ox . The 
results the analysis of the entire sample of 64 QSOs 
using the EM method available in ASURV (Isobe, 
Feigelson, & Nelson 1986) are shown below: 

a ox = (-1.662 ± .028) + (.042 ± .027) (z - 4) 

+ (.071±.019)(M B + 27) (1) 

= (-1.737 ± .023) + (.031 ± .029) (z - 4) (2) 

= (-1.648 ± .049) + (.068 ± .019)(M B + 27) (3) 

The coefficients for the regression analysis with 
the Buckley-James (B-J) and Schmitt's methods 
(if applicable) that are included in ASURV are 
consistent with those from the EM method shown 
above. 

From either of Eqs. 1 & 3, we find a 3.5cx de- 
pendence of a ox to Mb- The slope coefficeint can 
be translated into L x oc Lf pt with (3 — 0.56 ± 0.12 
(from Eq. 3), that is consistent with the recent 
extensive analysis including lower redshift QSOs 
by Strateva et al. (2005). 
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Fig. 2. — Grouping of QSOs from our sample and 
z > 4 samples from literature and (a ox ) of each 
group plotted against z and Mb- The QSOs are 
grouped according to the location in the z-Mb 
plane as plotted in panel (b). The groups are dis- 
criminated by different symbols as labeled. Radio- 
detected RLQs arc indicated by a large cross. The 
(a ox ) value of each group is plotted against z and 
Mb in panels (a) and (c) respectively. Error bars 
of (a ox ) are the standard deviation of the mean 
and those of z and Mb are the range of these val- 
ues in the group respectively. The data points for 
the groups are plotted with larger symbols than 
those for the individual QSOs. 
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Fig. 3. — The number densities of soft X-ray se- 
lected QSOs from samples used by SXLF3 (blue 
hexagons) and Silverman et al. (2005) (green 
squares) are plotted with those of optically se- 
lected QSOs (SSG95; Fan et al. 2001) (triangles), 
converted to our default cosmology. The X-ray lu- 
minosity cut (log L x > 44.95) was chosen to give 
the same number density as the optical curve at 
z ~ 2.9. The hatched area is a 90% confidence 
range of expected evolution of soft X-ray selected 
QSOs from the constraint on the systematic a ox 
change with redshift as determined from the com- 
parison between groups A & D. For reference, the 
density curves of QSOs from SXLF3 and Silver- 
man et al. (2005) have been overplotted for the 
original paper's representative luminosity ranges 
as labeled and plotted in open symbols. The dot- 
ted data points from SXLF3 show the absolute 
maximum number densities where all the optically 
faint unidentified X-ray sources were at the center 
of the redshift bin. 



Our primary goal of this study is to investi- 
gate the redshift dependence of a ox between z ~ 3 
and z > 4 to constrain the difference of number 
density behaviors of optically and X-ray selected 
QSOs at Mb ~ —26.5. In order to achieve this 
goal, we make more careful comparison of groups 
A (our sample) and D, which are well-separated 
in the redshift range ((z) of 2.93 and 4.79 respec- 
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tively) and occupy approximately the same lumi- 
nosity range (Fig. 2(b)). Since the difference in 
the mean Mb values for these two groups are small 
(= 0.16), the systematic difference of a ox between 
these groups due to optical luminosity dependence 
is negligible (~ 0.01 in (a ox )). The observed (a ox ) 
values of these groups are essentially the same. 
We have made Monte-Carlo simulations to find 
the probability distribution of Ada^ox), which 
is the difference between groups D and A to con- 
strain the systematic change in a ox with rcdshift 
at z > 2.7 at this luminosity range. In order to 
make a parent distribution of a ox from which the 
simulations draw objects, we used the all QSOs in 
the current sample except for upper limits. Also 
a ox of QSOs in each group have been shifted so 
that all groups have the same mean value. The 
standard deviation of these 58 a ox values is con- 
sistent with those of groups A and D respectively. 
In each run, we randomly took 6 QSOs to repre- 
sent group A and 11 to represent D respectively. 
The mean of each of the 6 and 11 random a ox 
values has been calculated and the distribution of 
the difference of these means (Ada(o! ox )) for 2000 
simulations have been investigated. 

As a result, the range where 90% of the simu- 
lations fall in were Ada(o:o X )da| < 0.13. This 
limits the systematic difference in log(L x /L pt) 
between z ~ 3 (group A) and z > 4 (group D) 
of ±0.34. For the slope of the luminosity func- 
tion (LF) of 7 X = 2, this corresponds to a num- 
ber density difference of |Alog(p pt/Px)| < 0.68, 
where p x and p opt arc the number densities of 
QSOs above optical luminosity and X-ray lumi- 
nosities respectively. Figure 3 shows the evolution 
of comoving number density of soft X-ray selected 
QSOs (SXLF3) as a function of redshift plotted 
with those of optical QSOs from SSG95 and Fan 
ct al. (2001) at z > 2.7 (their original values are 
converted to our default cosmology and a ). The 
luminosity limit for the soft X-ray QSOs is set at 
log L x > 44.95, where L K is the observed frame 
0.5-2 keV in erg s _1 ) corresponding to approxi- 
mately the 2 — 8 keV luminosity at the rest frame. 
Note that this is a rest-frame 0.5-2 keV luminos- 
ity, under the assumtion of T = 2.0 power-law 
spectrum, which is representative of X-ray selected 
type 1 AGNs and SXLF3 treats the luminosities as 
such. We also show the results from (Silverman et 
al. 2005)'s sample for the same observed 0.5-2 keV 



luminosity cut. This luminosity was selected such 
that the space density became equal to that of op- 
tical QSOs at z ~ 2.7 in thick lines and filled sym- 
bols. For reference, we show the number density 
curves from SXLF3 and (Silverman et al. 2005) 
(open symbols with thin lines), in lower luminos- 
ity ranges (see labels), which are representative 
of the respective samples. These show declines in 
z > 3. The large error bars in the log L x > 44.95 
data at z > 2.7 show that the QSOs in this regime 
is still underrepresented by current X-ray surveys. 
The limits of the number density curve of X-ray 
selected sample are shown in the shaded area in 
Fig. 3, corresponding to |A(a ox )DA| < 0.13. 

A limitation of the above investigation is the 
uncertainties in the effects of the sample selec- 
tion and variability. The X-ray and optical lumi- 
nosities have been measured in different epochs. 
Thus the variability of AGNs have a net effect 
of increasing the variance of the a ox distribution. 
Our underlying assumption is that the variabil- 
ity does not cause a net systematic difference in 
its effect on the mean a ox between sample A and 
sample D,separeted in redshift, but not in lumi- 
nosity. Both are optically selected samples and 
(a ox ) should be biased towards larger (more opti- 
cally luminous) values than the "true" (a ox ) (he. 
a ox of time-averaged mean optical and X-ray lumi- 
nosities) , because optical selection is more likely to 
pick up the AGN when it is more optically-bright, 
while the X-ray followup of the same object typi- 
cally gives average X-ray luminosity of the source. 
As long as both are selected in the optical and fol- 
lowed up by X-ray, that the slope of the LFs are 
the same at both redshifts, and that there is no 
systematic difference in the variability amplitudes 
of AGNs with redshift, the effect of this "variabil- 
ity" bias should be the equal between sample A 
and sample D. Thus we do not expect that the 
variabilty bias plays a major role in our analysis 
on the rcdshift dependence of (a ox ). 

From a combination of our sample at z ~ 3 
and z > 4 QSOs observed by Chandra, we have 
confirmed the apparent dependence of a ox on op- 
tical luminosity. This result, however, should be 
treated with caution, because our sample is op- 
tically selected and subject to the variability ef- 
fect of preferentially picking up optically brighter 
phase as described above. An X-ray selected sam- 
ple covering a much larger regime in z-L space 
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rather showed L x oc L opt (Hasinger 2005). This 
is, however subject to a similar effect working in 
the opposite sense. At this limited regime, on the 
other hand, there is a hint that our presently de- 
termined dependence might reflect the true behav- 
ior of the shift of a ox with QSO power. Expressing 
the LF as d$>/dlog L oc L~ 7 at the luminosities of 
interest, the soft X-ray LF has -fx = 2.2 ± 0.3 and 
(90% error) in 2.7 < z < 4.8 (from the same sam- 
ple as SXLF3). This can be compared with the 
optical LF of 7 opt — 1.6 ± 0.2 (la error) by Fan et 
al. (2001) or 7 opt = 1.87 by Schmidt, Schneider, & 
Gunn (1995). The trend that the optical LF has 
a flatter slope than the soft X-ray counterpart is 
consistent with the relation 7 opt = /?7 X within er- 
rors. This is also consistent with the comparison 
between X-ray and optically-selected AGN lumi- 
nosity functions by Ueda et al. (2003) (see their 
Fig. 20), where a conversion of their hard X-ray 
LF to the optical band assuming L x oc L op 7 t ° gave 
a good match to an observed optical QSO LF at 
high luminosities. However, our most recent com- 
parison of the SXLF (Hasinger, Miyajife Schmidt 
2005) (high luminosity end) and optical QSO LF 
by Croom et al. (2004) at z < 2.1 is more consis- 
tent with L x oc L opt , thus a more study is needed 
to investigate the relationship between direct com- 
parison of L x and L opt and the conversion between 
the X-ray and optically selected QSO LFs. 

5. Summary 

We have made Chandra ACIS-S observations of 
six QSOs at z ~ 3, which marks the peak of lu- 
minous QSO number density. These observatios 
fill a redshift gap in the X-ray coverage of lumi- 
nous QSOs. We found an average photon index 
of (r) = 1.9 ± 0.3 from the stacked spectrum and 
we also found (a ox ) = — 1.65 ± .05. The (a ox ) 
value is essentially the same as those at z > 4 in 
the similar UV luminosity range and thus we have 
found no systematic shift of X-ray to UV luminos- 
ity ratios with redshuft above z = 3. The density 
curves of Mb < —26.5 optically selected QSOs 
and log L x > 44.92 soft X-ray selected QSOs, giv- 
ing the same densities at z ~ 2.7, arc statistically 
consistent with each other within our limit of the 
systematic (a ox ) shift at z > 3. We note that 
this regime is still underrepresented by X-ray sur- 
veys. Large-area moderarely-deep X-ray surveys 
are needed to trace the rise of number density of 



the most luminous QSOs at z > 3 in X-rays. 
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Table 1 

Log of Chandra Observations and Sample Properties 



Name 


Obsid/Date 


Expo. 


z 










(ks) 








PC 0041+0215 


4150/2003 Sep 01 


9.2 


2.93 


19.5 


< 1.5 


PC 0947+5628 


4151/2003 Jan 25 


9.0 


2.91 


19.5 


< 1.0 


PC 1000+4751 


4152/2002 Dec 18 


13.9 


2.95 


20.0 


< 1.2 


PC 1015+4752 


4153/2003 Jan 01 


8.1 


2.92 


19.4 


< 0.9 


PC 1035+4747 


4154/2003 Mar 16 


10.0 


2.96 


19.6 


2.3 


PC 1447+4750 


4155/2003 Jul 28 


7.0 


2.93 


19.3 


< 0.9 



a Some authors prefer to use the notation "^4Bi45o(i+^)" • 
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Table 2 

X-ray Spectral Analysis Results and Derived Quantities 



Name 




AfG a 
JV H\20 


1 


cint b 
*x,14 


atz a 
JV H,20 


oCIS 


lOg Ll X 


Mb 


«ox 


PC 0041+0215 


2.93 


2.8 


2.4(1. 9;2. 9) 


1.4(1. 0;1. 8) 





36 


45.0 


-26.8 


-1.56 


PC 0947+5628 


2.91 


1.0 


1.5(0.7;2.3) 


.31(.17;.52) 





10 


44.4 


-26.8 


-1.78 


PC 1000+4751 


2.95 


0.9 


1.5(1.2;1.9) 


1.2(1.0;1.5) 





60 


45.0 


-26.3 


-1.48 


PC 1015+4752 


2.92 


1.0 


2.4(1. 8;3.0) 


.94(.65;1.2) 





24 


44.9 


-26.9 


-1.64 


PC 1035+4747 


2.96 


1.3 


0.6(-0.3;1.4) 


.17(.08;.34) 





9 


44.3 


-26.7 


-1.80 


PC 1447+4750 


2.93 


2.0 


2.0(1.4;2.6) 


1.0(.64;1.4) 





21 


44.9 


-27.0 


-1.65 


(6 QSOs) 


2.93 




1.9(1. 7;2.1) 




< 90 




160 







Note. — The 90% confidence ranges and upper limits for one interesting parameter are shown for free 
parameters. Spectral parameters without a confidence range were fixed during the fit. 

a In units of 10 20 cm- 2 . 

b X-ray flux before Galactic absorption in 0.5-2 keV (observer's frame), in units of 10~ 14 erg s _1 cm~ 2 . 
c Base-10 logarithm of X-ray luminosity at the rest-frame 2-10 keV band in units of erg s _1 . 



Table 3 

Additional Chandra Archival Data Analysis 



Name 


z 


Obsid/Datc 


Expo. 


JV H,20 


cint a b 


AB1450 


Rl 


M B 








(ks) 












PC 0910+5625 


4.04 


4821/2004 Mar 28 


23. 


2.9 


1.7(1.0;2.7) 


20.7 


< 1.5 


-26.2 


-1.67 


SDSS J235718.36+004350.3 


4.36 


4827/2003 Nov 26 


12. 


3.3 


5.0(3.3;7.2) 


20.2 


< 1.2 


-26.9 


-1.56 


SDSS J144428. 67-012344.1 


4.17 


4826/2004 Jan 8 


10. 


4.0 


1.8(0.7;3.4) 


19.8 


< 1.2 


-27.1 


-1.79 



a See notes to Table 2 for units. 
b 90% confidence ranges are shown. 



Table 4 
Comparison of a, 



group 




{*) 




(M B ) 


N Q 


(a ox ) a 


A 


2.9,3.0 


2.93 


-27.0,-26.3 


-26.75 


6 


-1.652±.046 


B 


3.5,4.6 


4.29 


-28.7,-27.7 


-28.35 


15 


-1.745±.042 


C 


3.5,4.6 


4.26 


-30.3,-28.7 


-29.22 


17 


-1.814±.035 


D 


4.0,5.3 


4.79 


-27.5,-26.0 


-26.91 


11 


-1.654±.031 


E 


4.6,6.4 


5.13 


-28.8,-27.5 


-27.98 


14 


-1.691±.047 



a Thc errors are the standard deviation of the mean. 
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